Robust optimal quantum gates for Josephson charge qubits 
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Quantum optimal control theory allows to design accurate quantum gates. We employ it to design 
high-fidelity two-bit gates for Josephson charge qubits in the presence of both leakage and noise. 
Our protocol considerably increases the fidelity of the gate and, more important, it is quite robust 
in the disruptive presence of 1// noise. The improvement in the gate performances discussed in this 
work (errors ~ 10~^ -7- 10"'* in realistic cases) allows to cross the fault tolerance threshold. 



One of the fundamental requirements of any proposed 
implementation of quantum information processing is the 
ability to perform reliably single- and two-qubit gates. 
In the last decade there has been an intense experimen- 
tal and theoretical activity to realize suitable schemes 
for quantum gates in a variety of physical systems as 
NMR, ion traps, cold atoms, solid state devices, just to 
mention a few [ij. Typically, as compared to single-bit 
gates, two-qubit gates are much more difficult to real- 
ize. The interaction between the qubits is more delicate 
to control while preserving coherence. Furthermore two- 
bit gates are more sensitive to imperfections, noise and, 
whenever present, leakage to non-computational states. 
It is therefore of crucial importance to find strategies to 
alleviate all these problems. A powerful tool to realize 
accurate gates is quantum optimal control 0], already 
applied for example to quantum computation with cold 
atoms in an optical lattice 3] . Aim of the present work is 
to apply optimal control to the realm of solid-state quan- 
tum computation, more specifically to qubits realized 
with superconducting nanocircuits. Josephson-junction 
qubits [1, IBI are considered among the most promising 
candidates for implementing quantum protocols in solid 
state devices. Due to their design flexibility, several dif- 
ferent versions of superconducting (charge, flux, phase) 
qubits have been theoretically proposed and experimen- 
tally realized in a series of beautiful experiments Q . Sev- 
eral schemes for oubit coupling have also been proposed 
(see the reviews [j, [1]). On the experimental side, cou- 
pled qubits have been realized in the charge 0, Q and 
in the phase 0] regimes where a CNOT and a ViSWAP 
gates have been implemented respectively. In the experi- 
ment of Steffen et al. Q the measured fidelity was of the 
order of 75% increasing up to 87% after accounting for 
measurement errors. Further improvements in the accu- 
racy rely on achieving larger decoherence times. In the 
experiment of Yamamoto et al. 0] a direct determination 
of the fidelity from the data was not possible, but it has 
been estimated to be ~ 80%. Advances in fabrication 
techniques will play a crucial role in achieving accurate 
quantum gates, however as the thresholds for fault toler- 



ant computation [10| are quite demanding, gate optimiza- 
tion is a powerful tool for a considerable improvement of 
their accuracy. A major open question is the resilience 
of optimized operations to imperfections affecting a real 
laboratory implementation, including: leakage to states 
outside the Hilbert subspace employed for logical encod- 
ing; inaccurate realization of the desired pulse shape; and 
classical noise in the system. 

In this Letter we apply optimal quantum control to su- 
perconducting charge qubits (that we choose for illustra- 
tion purposes). We analyze in detail the effect of noise 
and leakage, and we show that optimization keeps yield- 
ing a considerable improvement in gate fidelities even un- 
der such realistic conditions. In the context of supercon- 
ducting charge qubits, it has been proposed to couple the 
qubits via a capacitance 0, [U , an additional Josephson 
Junction (J J) or an inductance 13, The two-bit 
gate is realized by an appropriate choice of pulses in the 
gate potentials. For the two cases of capacitive and JJ 
coupling we construct the optimal pulse shapes thereby 
obtaining very high fidelities. For the case of capacitive 
coupling optimal control has been applied to sup ercon- 
ducting qubits for the first time by Sporl et al. [15|]. Here, 
we extend their results in two important aspects: First, 
we compare two different couplings in order to optimize 
the design. Second, we include the effect of 1// charge 
noise, believed to be the main source of decoherence in 
these systems [16|, ll7| , and show that the optimal gates 
are robust against it. We further show that gate accu- 
racy is maintained even under partially distorted pulse 
shapes. 

Coupled Josephson qubits Josephson charge qubits, 
sketched in Fig. [l] are defined in the regime in which 
the Josephson coupling is much smaller than the charg- 
ing energy. The single-qubit Hamiltonian (including also 
non-computational states) is defined as [1, Q 
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FIG. 1: a) A Cooper pair box can implement a charge qubit 
when tuned in the regime in which only two charge states are 
relevant. The box between the qubits represents the coupling 
which is specified below, b) An extra Josephson junction (left) 
or of a capacitance (right). 



where rii is the number of excess Cooper pairs on the 
i-th {i = 1,2) qubit, nf = CgVf^ l{1e) is the off- 
set charge controlled by the gate voltage Vg^^ (Cg is 
the gate capacitance), Ec is the charging energy and 
Ej is the Josephson coupling. By projecting onto the 
Hilbert space spanned by the states |0) , |1) {D = 2, 
D is the dimension of the Hilbert space) one recov- 
ers the charge qubit Hamiltonian. We want to include 
the effect of leakage to the charge states (in this case 
D > 2). Since we have Ej/Ec ^ 1, it is sufficient 
to add few other charge states. We included the charge 
states from |— 2) to |3)), i.e. D — 6. However in the range 
Ej /Ec ~ 5-^ 10 X 10^^ we verified that retaining the 
charge states | — 1) , |0) , |1) , |2) is sufficient. 

The coupling between the qubits (see Fig[T|3) can be 
either via a capacitor or a Josephson junction. In the 
case of capacitive coupling, Fig[T)3 (right), the interaction 
Hamiltonian reads 

Hf^ Ecc ^{ni - ng^i){n2 - rig, 2) ^12) {ni,n2\ (1) 
"1 ,"2 

where E^c is the charging energy associated to the 
Coulomb interaction between the qubits. If instead the 
coupling is via a Jospehson junction, FiglT|D (left), the 
coupling Hamiltonian is given by 

=^ ^(|ni)K + l)(ni + l|(n2K/i.c.) (2) 

where Ejj is the Josephson energy of the coupling junc- 
tion 0. 

Two- qubit gates The goal is to implement the univer- 
sal two-qubit gates Gjj and Gcc for the JJ and capacitive 



couplings respectively. They read 
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where we used the basis {|++) , |H — ) , | — h) , | )} for 

Gjj and the basis {|11) , |10) , |01) , |00)} for Gcc (|±) = 
(|0) ± |l))/\/2). Even under ideal operating condi- 
tions these gates cannot be implemented exactly [1, ■ 
As discussed in [l2j, Gjj can be approximately real- 
ized by tuning both qubits to degeneracy, fixing all the 
Josephson couplings to be equal in magnitude and turn- 
ing on the interaction for a time tjj ~ 0.97 27r/ Ejj. 
For Gcc we choose the same parameters of the ex- 
periment [| {Ej/E^^^ ~ 0.0777,Ej/eI?^ ~ 0.0610, 
Ecc/E^^^ — 0.1653). The time needed for the gate is 
Tec - 1-18 tt/e'-j''. Defining {a = JJ,cc) as the 
time evolution operator associated to the full Hamilto- 
nian Til + Ti.2 + Tif, a figure of merit to quantify the 
accuracy of a quantum gate is the error defined as 

e = l-Ti:{GiU?) (4) 

The Uis the time evolution operator projected onto the 
computational states (the fidelity of the operation should 
be tested only on the computational basis (|ni,n2) = 
|00),|01),|10),|11)). The fidelity is defined as = 1-e. 
In the following we determine optimized fidelities, up to 
a global phase for the gates in Eq. ([3]) when implemented 
with charge qubits. In order to be as close as possible to 
the experimental situation, we search for optimal pulses 
with the constraint that after the gates the two qubits 
are in their idle points. 

Optimized quantum gates Quantum control tech- 
niques described in 0, y| allow to minimize the error 
e defined in Eq. One assumes that the Hamiltonian 
is controlled by a set of external parameters which can 
be varied in time. The goal is to find the time depen- 
dence of the parameters which minimizes e. To illustrate 
it in a little more detail, let us imagine a system gov- 
erned by the time dependent Hamiltonian TClgit)], where 
g{t) is the control parameter. The goal of a quantum 
optimal control algorithm in general is to reach, in a cer- 
tain time r, a desired target state \iPt) with high fidelity. 
The algorithm employed here, due to Krotov 0, works 
as follows: (i) an initial guess go{t) is chosen for the con- 
trol parameter; (ii) the initial state |V'o) is evolved in 
time according to the dynamics dictated by Ti.[g(t)] un- 
til time r: lipgir)) — Ur{g)\'ipo)'j (iii) an auxiliary state 
IXg(T)) = |i/)t)(V't !'/'£(''')) is defined, which can be inter- 
preted as the part of \ipg{T)) that has reached the target 
H^t)', the auxiliary state is evolved backwards in time un- 
til t = 0; (iv) \xg{t)) and |V'g(*)) are propagated again 
forward in time, while the control parameter is updated 
with the rule git) ^ g{t) + Im [{xgit)\dgH\^git))] /X{t). 
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FIG. 2: (color online) Error Emm as a function of leakage 
for two-qubit gates: a) for Josephson-j unction coupling (in- 
cluding a residual capacitive coupling with Ecc/ Ejj = 0.05), 
with and without optimization (lower and upper curve re- 
spectively); b) for capacitive coupling, with optimization, as 
a function of the ratio E^p /Ecc (which we use here since the 
two charging energies in the experiment are different). 



The experimental value E'^j^ jE^ 



0.47 Ref. 



is marked. 



The weight function \(t) constrains the initial and fi- 
nal values of the control parameter; (v) steps (iii) and 
(iv) are repeated until the desired value of the fidelity 
is obtained. The same procedure can be followed also 
when the Hamiltonian contains more than one parame- 
ter. After a sufficient number of iterations, the algorithm 
converges and reaches asymptotically a minimum £m,in- 
In the present case, we consider E'f\ng^i, Ea as control 
parameters (Josephson couplings can be tuned by means 
of an applied magnetic fiux), and we look for optimal 
pulse shapes to improve the fidelity JF„. Although in 
principle one may consider all the different couplings in- 
dependently, this is impractical for an experimental point 
of view. In the case of JJ coupling we keep the gate 
voltage fixed and consider the same time dependence for 
all the Josephson couplings. This type of control can 
be achieved by applying a uniform time-dependent mag- 
netic field. In the case of capacitive coupling we allow for 
time-dependent gates but keep the Josephson couplings 
fixed. Relaxing these constrains will certainly lead to a 
further optimization of the fidelity at the cost, however, 
of a more complex external control. The important point 
is that already at the level discussed in this work the im- 
provement in the gate performances allows for crossing 
the fault tolerance threshold (lo| . 

The presence of leakage may be disruptive for two-bit 
gates in Josephson charge qubits fl9^| . Optimization how- 
ever fully compensates for leakage in both of the schemes 
depicted in Fig. [1] In the case of JJ coupling, Fig. [2] 
(left panel), we have only one control parameter, the 



Josephson coupling energy {E''j\t) = Ey'{t) = Ejj(t)). 
The non-o ptim ized gate (white circles) is realized as de- 
scribed in [12l | while the optimized curve, for the qubits 
of Ref. 01 {Ej/Ec ~ 3 X lO'^), gives an error of the or- 
der of 10"'*. This error is not appreciably influenced by 
the choice of the initial pulse, but rather it is physically 
determined by the constraints imposed on the pulse itself 




FIG. 3: Optimized gate error £„iin under noise having power 
spectrum Sng (i^) = A/lo as a, function of noise strength A: a) 
for Josephson coupling with Ej/Ec ~ 0.05, with and without 
optimization (lower and upper curves respectively); b) for ca- 
pacitive coupling with E^j ^ /Ecc ~ 0.47. Typical experimental 
value are around A ~ 10 



- for instance, requiring it to start and end at an optimal 
working point away from degeneracy, as we do here. In 
both cases we include leakage and the small effect of a 
finite charging energy Ecc- In the case of ca pac itive cou- 
pling, we build on the results obtained in [l5| and use 
their pulse sequence as the inial guess. Thus we present 
here only the optimized gate. Our results are shown in 
Fig. [2] (right panel) . In this setup, which coincides with 
that of the experiment of Ref. [81] the coupling Ecc can- 

should be chosen properly in 

Consistently, if e'^j^ jEcc 
Tec should be divided by 



not be changed. The values of the parameters i?j ■* / E)^ 
Ef/Ecc (i = 1,2) and 
order to realize the gate Gcc 
is changed by a given factor 



lEc 



the same factor. For the experimental value of Ej 

the error is £min — 10~^. Note that increasing E^f' jEcc 
results in a faster gate, thus reducing the effect of deco- 
herence. Here, in the best case, we can reduce the gate 
time to r 30ps, while keeping the fidelity constant. 
An important question to be addressed is to what extent 
our optimized gates are robust against noise. For this 
reason we check how stable the fidelity (optimized in the 
absence of noise) is, when the environment is taken into 
account [20|. The most important source of decoherence 
in charge qubits is 1// charge noise [2l[ . Although its un- 
derstanding is far from complete, 1// noise is believed to 
originate from two-level fluctuators present in the sub- 
strate and/or in the insulating barrier. Several theo- 
retical works have recently studied the relation between 
1// noise and decoherence in charge qubits (see 1^ 17 1 
and references therein). Here we follow the approach of 
Ref. [TtI I and model the environment as a superposition 
of bistable classical fluctuators resulting in an additional 
random contribution Sn^g\t) to the gate charge. A dis- 
tribution of switch rates 7 behaving as P{"f) oc I/7 in 
a range [jmin ■ "Jmax] results in a noise power spectrum 
Sn.iuj) = {Sn';}\t)Sn'^g\0))^ - uj-\ Following ^ we 
chose the switching rates such that the typical frequency 
of the gates is centered in between the two orders of mag- 
nitude over which the 1// noise extends ( We checked the 
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FIG. 4; Gate error Smin as a function of the pulse spectral 
cutoff uic for Josephson coupling with Ej/Ec = 1/20 (left) 
and for capacitive couphng with E^p /Ecc ~ 0.47 (right). In- 
sets: Corresponding optimal pulses. 



stability of our results with the choice of 7mm and ^max, 
data not shown) . We considered up to one thousand inde- 
pendent fl-uctuators coupled weakly to the qubits and we 
assumed that the charge noise on the two separate qubits 
is uncorrelated. The results of our analysis are reported 
in Fig. [3) regardless of the coupling scheme, the fidelity 
turns out to be quite robust against noise. Moreover, the 
error rates remain orders of magnitude better than with- 
out application of the quantum control algorithm, even 
under significant noise strengths, up to A ~ 10~^-^ 10"'^. 
We checked these results also with different kinds of noise 
(white noise, homogeneous frequency broadening in the 
control pulses) and we found similar conclusions (see 
also [iBl)- We finally investigated the dependence of the 
gate error on the experimentally unavoidable inaccura- 
cies of the pulse shapes. To this end we applied a filter 
to suppress the contribution of harmonics above a cutoff 
ijjc in the shape of the optimal pulses. In Fig. |4] we show 
the dependence of the error on the number of frequen- 
cies that compose the optimized pulses. In both cases 
the most important corrections are those at lower fre- 
quencies, as already pointed out in [l^. This explains 
the robustness of both optimized gates against noise pro- 
cesses: the fidelity is just marginally influenced by new 
frequencies introduced by the noise. Altough the realiza- 
tion of (nearly) optimal pulses is demanding, it definitely 
leads to accurate gate operation. One can then imagine 
to realize the two-bit gates in a longer time, in which 
case the shape of the pulse should be easier to realize. 
On the other hand, if the gate is too slow decoherence 
becomes relevant. It is then important to find an optimal 
gate time for which these two competing effects are mini- 
mized. We believe that this may be an avenue to realizing 
high-fidelity computations with Josephson nanocircuits. 

We thank the authors of Ref.[l3| for sharing with us 
information about the optimal pulses obtained in their 
work. We ackowledge support by EC-IHP (QOQIP) and 
EC-FET/QIPC (ACQP, EUROSQIP and SCALA), by 
the NSF through a grant to ITAMP at the Harvard- 
Smithsonian Center for Astrophysics and by Centro di 
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